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Abstract. In this article, we consider quantum crystals with defects in the 
reduced Hartree-Fock framework. The nuclei are supposed to be classical par- 
ticles arranged around a reference periodic configuration. The perturbation is 
assumed to be small in amplitude, but need not be localized in a specific re- 
gion of space or have any spatial invariance. Assuming Yukawa interactions, we 
prove the existence of an electronic ground state, solution of the self-consistent 
field equation. Next, by studying precisely the decay properties of this solution 
for local defects, we are able to expand the density of states of the nonlinear 
Hamiltonian of a system with a random perturbation of Anderson-Bernoulli 
type, in the limit of low concentration of defects. One important step in the 
proof of our results is the analysis of the dielectric response of the crystal to 
an effective charge perturbation. 

1. Introduction 

In solid state physics and materials science, the presence of defects in materi- 
als induces many interesting properties, such as Anderson localization and leads 
to many applications such as doped semi-conductors The mathematical modeling 
and the numerical simulation of the electronic structure of these materials is a 
challenging task, as we are in the presence of infinitely many interacting particles. 

The purpose of this paper is to construct the state of the quantum electrons of 
a mean-field crystal, in which the nuclei are classical particles arranged around a 
reference periodic configuration. We work with the assumption that the nuclear 
distribution is close to a chosen periodic arrangement locally, but the perturbation 
need not be localized in a specific region of space and it also need not have any 
spatial invariance. To our knowledge, this is the first result of this kind for Hartree- 
Fock type models for quantum crystals, with short-range interactions. By studying 
precisely the behavior of our solution, we are then able to expand the density of 
states of the Hamiltonian of the system in the presence of a random perturbation of 
Anderson-Bernoulli type, in the limit of low concentration of defects, that is when 
the Bernoulli parameter p tends to zero. The state of the random crystal and the 
mean- field Hamiltonian were recently constructed in [8]. Our small-p expansion is 
the nonlinear equivalent of a previous result by Klopp [19J in the linear case. 

The mean-field model we consider in this paper is the reduced Hartree-Fock 
model [3^, also called the Hartree model in the physics literature. It is obtained 
from the generalized Hartree-Fock model f25l by removing the exchange term. As 
the Coulomb interaction is long-range, it is a difficult mathematical question to de- 
scribe infinite systems interacting through the Coulomb potential. In the following. 
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we assume that all the particles interact through Yukawa potential of parameter 
TO > 0. In fact, we can assume any reasonable short-range potential, but we con- 
centrate on the Yukawa interaction in dimension d e {1,2,3} for simplicity. We 
consider systems composed of infinitely many classical nuclei distributed over the 
whole space and infinitely many electrons. 

We start by recalling the definition of the reduced Hartree-Fock (rHF) model 
for a finite system composed of a set of nuclei having a density of charge i^^^c and 
N electrons. The electrons are described by the A^-body wave-function (called a 
Slater determinant) 

ipixi,--- ,xn) = -^det{(pj{xi)), 

where the functions tfi E L^(R'^) satisfy {Lpi,(pj) — 6ij. The rHF equations then 
read 

H = -^A + V Vl<i<7V, (1) 

-AV + m^V = l^'^-i] - J^nuc) 

where pjp{x) = J2iLi ^^d Ai, • • • , Ajv are the smallest N eigenvalues of the 

operator H, assuming that A^r < Ajv+i. Here, |S'''~^| is the Lebesgue measure of 
the unit sphere S"^-^ {\S°\ = 2, \S^\ = 2tt, \S^\ = An). The existence of a solution 
of ([1]) is due to Lieb and Simon |26| . 

In order to describe infinite systems, it is more convenient to reformulate the 
rHF problem in terms of the one-particle density matrix formalism |24| . In this 
formalism, the state of the electrons is described by the orthogonal projector 7 — 
Si^i \ fi){'Pi\ rank N and the equations ([1]) can be recast as 



(2) 

(^7 ^nuc) : 

where formally p-y{x) — j{x,x) and the Fermi level ep is any real number in the 
gap [Aat, Aat+i). 

For infinite systems, the rHF equation is still given by ([2]), but 7 is now an infinite 
rank operator as there are infinitely many electrons in the system. The operator 
7 needs to be locally trace class for the electronic density Pj to be well-defined in 

The rHF equation ^ was solved for periodic nuclear densities 




^pcr 

ken 



by Catto, Le Bris and Lions in [10], and periodic nuclear densities with local per- 
turbations 

ri{- - k) + V 



were studied by Cances, Deleurence and Lewin in |7j. We have denoted by TZ the 
underlying discrete periodic lattice. The corresponding Hamiltonians are denoted 
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by iJpcr and iJ^. Stochastic distributions, 



nuc 



keiz ken 



for instance, were treated in |8]. 

Our present work follows on from [ZllSllH]. We are going to solve the equation ([2]) 
in the particular case where 



where v-pcr is a periodic nuclear distribution so that the corresponding background 
crystal is an insulator (the mean- field Hamiltonian iJpor has a gap around ep), and 
V G L^nif i^"^) is a small enough arbitrary perturbation of the background crystal. 
The perturbation v needs to be small in amplitude locally, but must not be local 
or have any spatial invariance. 

The rHF model is an approximation of the iV-body Schrodinger model, for which 
there is no well-defined formulation for infinite systems so far. The only available 
result is the existence of the thermodynamic limit of the energy: the energy per 
unit volume of the system confined to a box, with suitable boundary conditions, 
converges when the size of the box grows to infinity. The first theorem of this form 
for Coulomb interacting systems is due to Lieb and Lebowitz in [22\. In this latter 
work, nuclei are considered as quantum particle and rotational invariance plays a 
crucial role. For quantum systems in which the nuclei are classical particles, the 
thermodynamic limit was proved for perfect crystals by Fefferman [12J (a recent 
proof has been proposed in [17 ) and for stationary stochastic systems by Blanc 
and Lewin [4 . Similar results for Yukawa interacting systems are simpler than for 
the Coulomb case and follow from the work of Ruelle and Fisher [13J for perfect 
crystals and Veniaminov [32J for stationary stochastic systems. Unfortunately, very 
little is known about the limiting quantum state in both cases. 

For (orbital-free) Thomas- Fermi like theories, the periodic model was studied 
in [261 [5], the case of crystals with local defects was studied in [5J and stochastic 
systems were investigated in [3]- To the best of our knowledge, the only works 
dealing with systems with arbitrary distributed nuclei are [9, 2J for Thomas- Fermi 
type models. 

As mentioned before, our work is the first one to consider this kind of sys- 
tems in the framework of Hartree-Fock type models. Our results concern small 
perturbations of perfect crystals interacting through short-range Yukawa potential. 
Extending these results to more general geometries and for the long-range Coulomb 
interaction are important questions that we hope to address in the future. 

After having found solutions of ^ for any (small enough) u E -^unif (^'^)^ '^^ 
study the properties of this solution for local perturbations ly. This enables us to 
investigate small random perturbations of perfect crystals. Precisely, we consider 
nuclear distributions 



nuc 



= I'pcr + 



(3) 




where {qk)ken are i-i-d. Bernoulli variables of parameter p and x is a compactly 
supported function which is small enough in L'^{M.'^). We are interested in the 
properties of the system in the limit of low concentration of defects, that is when 
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the parameter p goes to zero. We prove that the density of states of the mean- 
field Hamiltonian Hp = — + V^, which describes the collective behavior of the 
electrons, admits an expansion of the form 

,7 

np = n^ + Y.^'lP^ +0{p'+^). (4) 
j=i 

Here, no is the density of states of the unperturbed Hamiltonian H-pcr — — ^A + Vpcr 
and is a function of the spectral shift function for the pair of operators Hp^t and 
H~^, the latter being the mean- field Hamiltonian of the system with only one local 
defect constructed in ffj. We give in Theorem 12.71 a precise meaning of 0(p-'+^). 
In |19j . Klopp considers the empirical linear Anderson-Bernoulli model 

H=-]^A + Vq + V with V{LU,x)^^qk{Lu)r]{x-k), 

ken 

where Vq is a linear periodic potential and 77 an exponentially decaying potential. 
He proves that the density of states of the Hamiltonian H admits an asymptotic 
expansion similar to Q. The case where V{ijJ,x) is distributed following a Poisson 
law instead of Bernoulli is dealt with in [50]. Our proof of (HJ follows the same lines 
as the one of Klopp. The main difficulty here is to understand the decay properties 
of the mean- field potential V solution of the self-consistent equations ([2]). For 
this reason, we dedicate an important part of this paper to the study of these 
decay properties. In Theorem 12.31 below, we show that for a compactly supported 
perturbation the difference V — V-pcr decays faster than any polynomial far from 
the support of the perturbation v. Moreover, we show that the potential generated 
by two defects that are far enough is close to the sum of the potentials generated 
by each defect alone. 

The article is organized as follow. In Section[21 we present the main results of the 
paper. We start by recalling the reduced Hartree-Fock model for perfect crystals 
and perfect crystals with local defects in Section [2T] In Section [2T2l we state the 
existence of solutions to the self-consistent equations ([2]) for Vmic given by ([3]) . We 
also explain that our solution is in some sense the minimizer of the energy of the 
system. We also prove a thermodynamic limit, namely, the ground state of the 
system with the perturbation v confined to a box converges, when the size of the 
box goes to infinity, to the ground state of the system with the perturbation v. 
In Section [5751 we prove decay estimates for the mean- field density and potential. 
In Section 12. 4[ we present the expansion of the density of states of the mean-field 
Hamiltonian. The proofs of all these results are provided in Sections 01 [H IH] and [T] 
In Section [31 we study the dielectric response of a perfect crystal to a variation of 
the effective charge distribution, which plays a key role in this paper. 

Acknowledgement. I thoroughly thank Eric Cances and Mathieu Lewin for their 
precious help and advices. The research leading to these results has received funding 
from the European Research Council under the European Community's Seventh 
Framework Programme (FP7/2007-2013 Grant Agreement MNIQS no. 258023). 

2. Statement of the main results 

2.1. The rHF model for crystals with and without local defects. In defect- 
free materials, the nuclei and electrons are arranged according to a discrete periodic 
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lattice TZ of M"^, in the sense that both the nuclear density ly^uc = i^pa- and the 
electronic density are 72.-periodic functions. For simplicity, we take 7?. = Z'' in the 
following. The reduced Hartree-Fock model for perfect crystals has been rigorously 
derived from the reduced Hartree-Fock model for finite molecular systems by means 
of thermodynamic limit procedure in |10|[7] in the case of Coulomb interaction. The 
same results for Yukawa interaction are obtained with similar arguments. The self- 
consistent equation ^ then reads 

70 = 1 (i/por < ep) 

^^pcr--^A + Fper (5) 
-AVpcr + mVpcr = \S'^~^\ {pj, - Z/por) ■ 

It has been proved in jlOl [7] that ([5]) admits a unique solution which is the unique 
niinimizer of the periodic rHF energy functional. 

Most of our results below hold only for insulators (or semi-conductors). We 
therefore make the assumption that 

-ffpor has a spectral gap around ep. (6) 

The rHF model for crystals with local defects was introduced and studied in [7J. 
A solution of the rHF equation ^ is constructed using a variational method. One 
advantage of this method is that there is no need to assume that the perturbation 

is small in amplitude. The idea is to find a minimizer of the infinite energy of 
the system by minimizing the energy difference between the perturbed state and 
the perfect crystal. The ground state density matrix can thus be decomposed as 

7 70 + Qu, (7) 
where Q,y is a minimizer of the energy functional 

S'iQ) = Tr^„ ((i/pe,. - eF)Q) + \d^{pq -v,pq- v) (8) 
on the convex set 

/C - {q* - Q, -70 < < 1 - 70, (-A + 1)^ Q e ©2(i'(M'')), 



(-A + 1)^ Q±± (-A + 1)^ e 6i(i2(K'^))} , 



(9) 



where Q++ = (1 - 7o)g(l - 70), = 7oQ7o and Tr^„(^) = Tr (A++ + A—). 

We use the notation (3p to denote the p"^ Schatten class. In particular 62 is the set 
of Hilbert-Schmidt operators. The second term of (jS]) accounts for the interaction 
energy and is defined for any charge densities f,g E i/^^(M'') by 



Drnif,g) = \S^-'\ f M^clp^ f f f{x)Y^{x-y)g{y)dxdy, 
JK<i IpI + jRd jRd 

where f[p) = (27r) Jjjd / {x) e~^P'^dx is the Fourier transform of /. The Yukawa 
kernel Y„i, the inverse Fourier transform of |S'''~^| (|p|^ + rri^)^^ , is given by 

Y^{x) = { Koim\x\) ifd = 2, 
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where Kq (r ) — e^ '' * dt is the modified Bessel function of the second type [23] . 
It has been proved in |7] that the energy functional (jS]) is convex and that all its 
minimizers share the same density p^. These minimizers are of the form 

j^l{H <ep)+S 

H = -^A + V (10) 

~AV + = I 5^-1 1 {p^ - I/per - 1^), 

where < S < 1 {H = e^). If z/ is small enough in the i/~^-norm, then S = 0. 

One of the purposes of this article is to find decay estimates of the potential V 
solution of (ITU)) that are necessary in the study of the Anderson-Bernoulli random 
perturbations of crystals. 

2.2. Existence of ground states. In this section, we state our results concerning 
the electronic state of a perturbed crystal. The host crystal is characterized by a 
periodic nuclear density Vpa- G L'^nH (^'') ^^^h that the gap assumption ^ holds. 
The perturbation is given by a distribution i' <E iu„if (R''). The total nuclear 
distribution is then 

In Theorem 12.11 below, we show that if is small enough in the L'^^^^f -norm, then 
the rHF equation ^ admits a solution 7. This solution is unique in a neighborhood 
of 7o. The proof consists in formulating the problem in terms of the density and 
using a fixed point technique, in the spirit of |15| . 

Theorem 2.1 (Existence of a ground state). There exists ac > and C > such 
that for any v € L^^jf (R'^) satisfying \\i^\\j^2 < ac, there is a unique solution 

7 £ ©ijoc(-^^(R'^)) to the self- consistent equation 
1 = t {H < ep) 

H = -]^A + V (11) 

-AV + m^V ^ \S'^-^\ (p-,-^-^pcr) 
satisfying 

IIP7-P7oIIl2_^,,. <CM^,^^^^^ . (12) 

We denote this solution by the response electronic density by — p^^ — p-^^ 
and the defect mean- field potential by ^ V — Vpa- ■ 

For a local defect v e L'^{R'^) n L^{R'^) such that ||i^||^2 ^ < ac, equation (fTT|) 
admits a unique solution which coincides with the ground state 7 solution of ([7]) 
constructed in [7]- Indeed, the solution 7^ given in Theorem l2.1l is a solution of the 
defect problem PH)) . Moreover, in the proof of Theorem 12.11 we prove that H has 
a gap around ep, thus necessarily (5 = in (|10p . As all the solutions of PU)) share 
the same density, (jTU]) (thus PT|) ) admits a unique solution. 

The ground state constructed in Theorem 12. II is in fact the unique minimizer of 
the "infinite" rHF energy functional. Indeed, following ideas of [16J, we can define 
the relative energy of the system with nuclear distribution Vnnc by subtracting the 
"infinite" energy of ji, from the "infinite" energy of a test state 7: 

^'fil) ■= Tr^.^ {{H - ep) (7- 7<>)) + ^D„, {p^ - p^^,p^ ~ p^J . 
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This energy is well-defined for states 7 such that 7 — 7,^ is finite rank and smooth 
enough for instance, but one can extend it to states in a set similar to /C in ([9|). 
The minimum of the energy 5^°' is attained for 7 = 7,^ = 1 {H < ep). Moreover, 
as H has a gap around ep, S]^^ is strictly convex and 7^ is its unique minimizer. 

In the following theorem, we show that if we confine the defect v to a box of 
finite size, then the ground state of the system defined by the theory of local defects 
presented in Section [2.11 converges, when the size of the box goes to infinity, to the 
ground state of the system with the defect v defined in Theorem 12.11 We denote 
by Ti = [-L/2,L/2f. 

Theorem 2.2 (Thermodynamic limit). There exists etc > such that for any 
V e -^unif (K'^) satisfying ||i^||^2 < Uc, the sequence (7i/irj^ )LeN\{o} converges in 
©ijoc(i^(K'')) to -f^ as 00. 

2.3. Decay estimates. In this section, we prove some decay estimates of the mean- 
field potential Vi, and the mean- field density , which will be particularly important 
to understand the system in the presence of rare perturbations in the next section. 

Theorem 12.31 below is crucial in the proof of Theorem 12.71 Indeed, we will need 
uniform decay estimates for compactly supported defects, with growing supports 
and uniform local norms. 

Theorem 2.3 (Decay rate of the mean- field potential and density). There exists 
ac,C' > and C > such that for any v £ L'^(U.'^) satisfying ||i^||^2 < etc we 
have for R. > 2 

where Cii{v) = {x G M'*, d (x, supp(i/)) < R] . 

Remark 2.4. Using the same techniques as in the proof of Theorem \2.3l we can 
prove (see [21j ) that there exists a,ac,C' > and C > such that for any v £ 
Lc(I^'^) satisfying ||i^||^2 < etc and \\i^\\fj-i < ct, we have for R> 2 

IIK|1^2(r.\Ch(.)) + IIp^IIl-^(m^\c«(.)) < Ce-^'(^°s«)^ 1,^1,^^^^^^ . (14) 

Estimate (jl4p gives a decay rate of the solution of the rHF equation for crystals with 
local defects, far from the support of the defect. In particular, it shows that G 
L^(W^). This decay is due to the short-range character of the Yukawa interaction. 
In the Coulomb case, it has been proved in [3] that for anisotropic materials, pi, ^ 
Li(R'*). 

The decay rate of Vy and p^, proved in Theorem 12.31 is faster than the decay of 
any polynomial, but is not exponential, which we think should be the optimal rate. 

Proposition 12.51 below is an important intermediary result in the proof of Theo- 
rem [2?2l It says that the mean- field density pi, and potential Vi, on a compact set 
depend mainly on the nuclear distribution in a neighborhood of this compact set. 

Proposition 2.5 (The mean-field potential and density depend locally on v). 

There exists ckf, > such that for any P > 2 there exists C > such that for 
any v G I^^^si {W^) satisfying ||i'||^2 < etc and any L > 1 we have 

C 

where vl = i^^Fl ■ 
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In the same way, we obtain the fohowing result which wih be very useful in the 
proof of Theorem 12.71 We prove that the potential generated by two defects that 
are far enough is close to the sum of the potentials generated by each defect alone 
in the sense of 

Proposition 2.6. There exists etc > such that for any (3 > 2, there exists 
C > such that for any h'1,1'2 G L'^{M.'^) satisfying 11^^111^2 , ||i^2||/^2 < etc o,nd 
R = d(supp(i^i), supp(i'2)) > 0, we have 




Proof. The proof is the same as the one of Proposition 12.51 with v = vi + V2 and 
L = 2R. □ 

2.4. Asymptotic expansion of the density of states. In this section, we use 
our previous results to study a particular case of random materials. In the so-called 
statistically homogeneous materials, the particles are randomly distributed over the 
space with a certain spatial invariance. More precisely, the nuclear distribution 
(thus the electronic density) is stationary in the sense 

^nuc 

where {Tk)k<=:i'i is an ergodic group action of on the probability set (see 
Figure [Ij . One famous example of such distributions is the Anderson model 



• • • *® (zf 

Perfect crystal Statistically homogeneous material 



Figure 1 . Example of a stationary nuclear distribution 

i/„uc(w, a;) = ^ '7fe(w) x{x - k), 

where, typically, x £ C^(R'^) and the g^'s are i.i.d. random variables. The reduced 
Hartree-Fock model for statistically homogeneous materials was introduced in [8J. 
The state of the electrons is described by a random self-adjoint operator {'y{^))i^^fi 
acting on L^(]R'') such that < 7(w) < 1 almost surely. The rHF equation is then 

(7(0;) = l{H{ij) < eF) + S{ij) 
H{oj) = -^A + V{uj, •) almost surely, (15) 

-Al/(c., •) + mV(c., •) = \S''-^\ (p^(^) - v{,o, ■)) 
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where < 6{uj) < almost surely. The solutions of (fT5|) turn out to be 

the minimizers of the energy functional 

where Ti{A)=E (Tr (Ir^lr)) and 

D„^{f,g)=E( f f f{x)Yra{x~y)g{y)dxdy 



Here, T = [-1/2, 1/2) denotes the semi-open unit cube. Thanks to the convexity 
of it has been proved in [8 that the minimizers of S.^^^^^ share the same 

density. Therefore, the Hamiltonian H solution of (|15l) is uniquely defined. 

In this paper, we are interested in the particular case of random perturbation of 
perfect crystals 

fnuc(w, x) = l'por(2:) + x) 

in the limit of low concentration of defects. We restrict our study to Anderson- 
Bernoulli type perturbations, that is, we suppose that at each site of Z'', there is a 
probability p to see a local defect independently of what is happening in the other 
sites. More precisely, we consider the probability space il = {0, 1}^ endowed with 

the measure P = (pSi + (1 — p)Sq)'^^ and the ergodic group action rfe(a;) = w.+fe. 
The defect distribution we consider is then given by 

i^p{uj,x) = ^ qk{uj)xix-k) 

where qk is the k^^ coordinates of w and % S LP'iW^) with supp(x) C F. The 
gfc's are i.i.d. Bernoulli variables of parameter p. If ||xIIl2 < Oic, then 5{u)) — 
almost surely and (|15p admits a unique solution. For almost every w, this solution 
coincides with the solution of ([TT]) constructed in Theorem 12.11 For convenience, 
we will from now on use the notation 

where we recall that e^? is the Fermi level. We introduce the mean- field Hamiltonian 
corresponding to the system with the defect Vp 

Hp = Hq + with (i^, x) = Y„i * (p^p - Vp) . 

As Vp is stationary with respect to the ergodic group {Tk)k&'^ ^'^d uniformly 
bounded in O x E'^, then by |271 Theorem 5.20], there exists a deterministic positive 
measure np{dx), the density of states of Hp, such that for any tp in the Schwartz 
space S(R) 

Lp{x)np{dx) = Tr {ip{Hp)) . 



For K C Z"^, we define the self-consistent operator corresponding to the system 
with the defects in K 

Hk = Hi} -^Vk, 

where 

Vk = Ym* {pK - vk), x{- - k) and pK = Pv^ ■ 

k^K 
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If \K\ < oo, we denote by £,k{x) the spectral shift function |33J for the pair of 
operators Hk and Hq. It is the tempered distribution in 5'(R) satisfying, for any 
if € S(R), 

Tr {^p{Hk) - f{Ha)) = / £,K{x)tf' {x) dx = - / ^'ii{x)ip{x) dx. 



In Theorem p. 71) below, we give the asymptotic expansion of the density of states 
Up in terms of powers of the Bernoulli parameter p. 

Theorem 2.7. For x G such that supp(x) C T and K C such that 

\K\ < oo, we define the tempered distribution px by 

' ' K'GK 

There exists Uc > Q such that if HxIIl^ — '^c, then 

(i) for i G {1,2}, pLj = ^ jfcz'', t^K is a well-defined convergent series in 

\K\=]. oeK 

S'(R). 



(ii) for J < 2, there exists Cj > 0, independent of x such that for any ip Cz S{ 

<Cj\\x\\l^ E ^a,0i^)p'+\ 



J 

j=i 



Q<(J + 3)(d + l) 

3< J + 4+(J + 2)d 



where uq is the density of states of the unperturbed Hamiltonian Hq and 

J^a,f}{^) = sup 



In Theorem 12.71 we only present the expansion of the density of states until the 
second order J — 2. The proof of the expansion up to any order J G N should 
follow the same lines and techniques used here. 

A result similar to Theorem 12 . 71 was obtained in [T9] in the linear case. Materials 
with low concentration of defects were studied by Le Bris and Anantharaman yy . 
in the framework of stochastic homogenization. 

The proof of Theorem 12.71 follows essentially the proof of Theorem 1.1]. It 
uses the decay of the potential related to each local defect. In [151 Theorem 1.1], 
the linear potential is assumed to decay exponentially. In our nonlinear model, the 
decay estimates established in Section 12.31 play a crucial role in the proof. 

The rest of the paper is devoted to the proofs of the results presented in this 
section. In the next section, we study the dielectric response of the crystal to an 
effective charge perturbation. The results of Section[3]will be used in later sections. 

3. Dielectric response for Yukawa interaction 

In this section, we study the dielectric response of the electronic ground state of 
a crystal to a small effective charge perturbation / e -^u^if (M'^). This means more 
precisely that we expand the formula 

Qf = t{HQ + f*Y„,<Q)-t (Ho < 0) 

in powers of / (for / small enough) and state important properties of the first 
order term. The higher order term will be dealt with later in Lemma 14.11 For 
Coulomb interactions and local perturbation / € L'^iW^) nCo(M''), where C^iW^) is 
the Coulomb space, this study has been carried out in |6] in dimension d = i. 
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The results of this section can be used in the Hnear model or the mean-field 
framework. In the reduced Hartree-Fock model we consider in this paper, the 
effective charge perturbation \s f = — u, where is the electronic density of 
the response of the crystal to the nuclear perturbation v defined in Theorem 12.11 
Expanding (formally) Qf in powers of / and using the resolvent formula leads to 
considering the following operator 

1/1 1 

ZITT Jc Z — Hq Z — Hq 

where C is a smooth curve in the complex plane enclosing the whole spectrum of Hq 
below (see Figure[2]). By the residue Theorem, the operator Qij does not depend 



a(i/o) 

, , , , 



Figure 2. Graphical representation of a contour C C C enclosing 

(7(Ho)n(-00,0]. 

on the particular curve C chosen as above. We recall that Vpcr is —A bounded with 
relative bound 0. Thus Hq is bounded below by the Rellich-Kato theorem |28[ 
Theorem X.12]. Theorem 13 . 1 1 below studies the properties of the dielectric response 
operator £ : / — > pg^ ^. and the operator (1 -I- £) , which will play an important 
role in the resolution of the self-consistent equation ([TT|) . In particular, it gives 
the functional spaces on which C and (1 + are well-defined for both local and 
extended charge densities. It also says that (1 + C)^^ is local in the sense that its 
off-diagonal components decay faster than any polynomial. We consider H~^{M.'^), 
endowed with the scalar product 

,r . 1 /• W)9iP) , 

{^-^r Jr-' \p\ + m? 
Theorem 3.1 (Properties of the dielectric response). We have 

(i) The operator 

f ^ -PQi.f^ 

is well-defined, bounded, non-negative and self-adjoint. Hence 1 + C is 
invertible and hicontinuous. 

(ii) The operator C is hounded from H-'^{W^) to L'^{W^) and 1/(1 + C) is a 
well-defined, hounded operator from i^(M'') into itself. 

(iii) The operator 

C: iS„if(K') ^ iLf(R') 

is well-defined and hounded. The operator 1 + C is invertihle on iy^if (M'') 
and its inverse is bounded. 
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(iv) There exist C > and C > such that for any j, fc e Z"* such that 
\k — j\ > 1, we have 

1 



r+i 



:1 



r+k 



(16) 



Proof. The proof consists in the fohowing 6 steps. In the whole paper C > and 
C" > are constants whose value might change from one hne to the other. 
Step 1. Proof of (i). The proof is similar to the one of [HI Proposition 2], with the 
Yukawa kernel y,„. instead of the Coulomb kernel. In the Yukawa case, H~^(M.'^) 
plays the role of the Coulomb space. The proof of ^ Proposition 2] can easily be 
adapted to our case. We skip the details for the sake of brevity. 
Step 2. Proof of (ii). Let / e H-^[W^). Then Y„, * / e L'^{m/) and 



\Y^*f\\ 



I 



L2 



jdp < c 



fip) 



\p\ + 



dp = C\\f\\i-, 



(17) 

Therefore, by |^ Proposition 1], Qij G IC, where /C has been defined in and 
Cf = -pQij G L'^(R'^). Arguing by duality, we have for any W G L'^{R''), 



Tr iQijW) = [ pQ.^^W. 
Besides, by the Kato-Seiler-Simon inequality [30l Theorem 4.1] for d < 3 



Vp > 2, ||/(-zV)5(x)||e^ < i2nr^ ||/||^, ||g|| 



LP 



and the fact that 



(z — Hq) ^ (1 — A) is uniformly bounded on the contour C, 



(18) 

(19) 
(20) 



we have 



^ Ym*f-^We&2{L^iR'')) 



z-Ho 



and 

|Tr (QijW)l 



1 

2i7r 



Tr 



z~Ho 



^ Yrn * f—^W I dz 



Z~Hq 



z~Hq 



<C\\Y^*fh, \\W\\^,. 

(21) 



The bound (I^Ul) follows from the following lemma. 

Lemma 3.2. Let W £ i^^^jj (R'^). Then there exists C > 0, depending only on the 
-^unif -iT-orm of W , such that for any z € C \ cr(— A + W), we have 

||(-A + 1)(-A + W~ z)'i|L < C- l^l""^ 

' ~ d(z,CT(-A + VF)) 

In particular, if K isacompactsetofC\(T{—A + W), then (—A + l){—A + W — z)~^ 
is uniformly bounded on A. 

The proof of Lemma [X^ is elementary, it can be read in pT]. In view of (|17p . H 
and dUl), it follows that 



{Cf)W 



<C\\f\\H-i\\W\\ 



L2 
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We deduce that 

<c||/||^_,. 

We now prove that (1 + Cy^ is bounded on L'^{M.'^). Let g e L'^{W^) and / e 
H-^iW^) such that (1 + £)/ = g. Then, f = g - Cf e L^iR'^). As 1/(1 + C) is 
bounded from H^^{M.'^) into itself, we have 

UWh-^ <C\\g\\H-^ <C\\g\\^.. 
Therefore, as C is continuous from H-^{R'') to L'^{R'^) , 

II/IIl^ = 11.9 - ^/IL. < Ml-^ + WCfh-^ < Ml-^ + C 11/11^-, < C \\g\y , 
which concludes the proof of (ii). 

Step 3. Proof of the first part of (in): C is well-defined and bounded on i^^^jj (R'^). 
First, we consider a bounded operator A S B{L'^{M.'^)) and prove that {z-Ho)-'^A{z- 
Ho)~^ is locally trace class. For x G L'^(R'^) and z G C, we have by ((20)) and the 
Kato-Simon-Seiler inequality (fT9l) that x{z — Ho)^^ A{z — Ho)^^x is trace class and 
that there exists C > independent of z £ C such that 



Tr X 



Ho 



-A- 



Hn 



< 



< 



1 



1 



-A- 



1 



z - 



Ho 



z-Ho' 



b2 



Ho 



<C\\A\\^\\x\\i.. 



It follows that the operator (z — Ho)^^A{z — Ho)^^ is locally trace class and that 
its density pz is in Ll^^{R'^). We now show that pz is in fact in i^^^jj. (R''). Let 
fc e Z'' and M be a non- negative function in L°°{r + k). It holds, taking x = V^j 
that 



PzU 



PzX 



Tr X 



1 



-A- 



1 



-X 



<C\\A\\^\\u\\^,. (22) 



z - Ho z - Ho 
By linearity, we deduce that pz £ L°°{R'^) and 

As all these estimates are uniform on the compact set C, the operator 

(2?7r)~^ {z — Ho)~^ A{z — Ho)~^ dz is locally trace class and its density p is in 

L^^jif {W^) and satisfies 

MlI^^^^ <C\\A\\^. (23) 

We now consider the case when A = y,„ * / is a potential generated by a charge 
density / G i^^if (R'^). The following Lemma gives the functional space Y,n * / 
belongs to when / G L^^jf (M'^). 



Lemma 3.3. Let / G i^^jf (K'^) and Y G if„, 



(24) 



for some 1 < p,q < oo. Then, the function Y * f is in i^^if {R'^) with 1 + 1/?' 
1/p + 1/q and there exists C > independent of f such that 

\\Y*fh. ^ <C\\f\\^. . 
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The proof of Lemma [331 is exactly the same than the one of [8j Lemma 3.1], we 
omit it here. As Ym satisfies ([M)) for p = 2, we have 



Therefore, by ([M]) 



and \\Y^*fh^<C\\f\\^. 



(25) 



<C||K„,*/||^„ <C 11/11 



which proves that C is weh-defined and bounded from L^^^^ 
concludes Step 3. 



l"^) into itself. This 



In the rest of the proof, we use a localization technique. We will thus need Lem- 
mas |331 and |33] below. Lemma 13.41 gives an estimate on the commutator between 
the dielectric response operator C and a localizing function in both L'^{R'^) and 
^unif i^'^)- Lemma [3.51 gives a decay rate of a real sequence satisfying a recursion 
relation that will be satisfied by the localized sequence. The proofs of Lemmas 13.41 
and 13.51 are postponed until the end of the proof of the proposition. 

Lemma 3.4. Let x be a smooth function in C^{M.'^) such that 0<x<l,X=l 
on B{0, 1) and x = outside B{0, 2). For any set / C Z'' and R> 1 we denote by 
Bi^R = Ufcg/ (S(0, R) + k) and by Xi,r{^) ~ X (d(a;, /) jB) . The family of functions 
{xim)r>i satisfy < xiM < 1; XiM = ^ on Bj^r, xi.r = outside Bi^2R and 



R\Vxi,rM+R^\^XiAx)\ <C fl.e.. 



(26) 



where C is independent of the set I . We denote by rji,R = 1 — Xi,R- Then, there 
exists C > and C" > such that for any I d 2,'^ and any f g L^(R'^j3 

hl,RYm * f -Yjn* {VlMf)\\H2 + WiVlM,^ /IIl2 



- R\ 



< 



c 

R 



e-^'^ll/IL. 



^Bi ,3r\Bj 



L2 



H-1 

(27) 



and for any f G Ll^,, [W') 



C'R I 



/II 



Bl .3r\Bj /2 



(28) 



Lemma 3.5. Let (a;ij)ij>o be a non-increasing family of real numbers such that for 
any R > 0, 

%Rla (29) 



XR < ^e'^"'^xo 
it 



R 



for given C > and C",a > 0. Then, there exists C > and C > such that for 
any R>2 

a;fl<Ce-^'('°s«)'a;o. (30) 



We now proceed with the proof of Theorem 13.11 We first prove (iv), then we 
prove that 1 + £ is invertible on L^^jj {W^). 



the whole paper, we use the convention f * gh = hf * g = h{f * g). 
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Step 4. Proof of (iv). We explain how to use Lemmas 13. 4l andl 3.5l to prove (|16|) . Let 
keZ''^ and for i? > 1, let ?7fl ~ '']{k},R a-iid _Bfl = -B{fe},_R as defined in Lemma p.4[) . 
Let g e L2(]^rf) a^j^(j denote by / = (1 + /:)"^ lr+fc5- For i? > f , we have 

riR (/ + Cf) = riRlr+kg = 0. 

Therefore 

(1 + C) r]Rf = f]Rf + C-qnJ = Cijuf - rjnCf = [£, -qn] /. 
Since f/(l + £) is bounded on L'^{W^), it follows that 
1 



<c\\[c,ini,]f\\ 



L2 



L2 



c 



<^e-^^ 11/11,. 



C_ 
R 



L2 



(31) 



where we have used Lemina l3.4l in the last step. Denoting by xji = ||1r<'\_b,;,/||,2j 
the estimate (PT|) leads to 

C C' R ^ 

XR< —e xq + —xii/4. 
K It 

Therefore, Lemma [3.51 gives that there exists C > and C" > such that for any 
R>2 

h^/IL. < X,,/, < Ce~^'(>°^«)^xo = ||/||,. < Ce-^'i^^^^)' hh^^r+k) 

where the last inequality follows from the fact that is bounded on L^(IR'*). 

Finally, as tr+j < V\k-j\/i-i/2, then 

^ ' < Ce-c^'(i°gl'=-^D' Mmr+k) < ^'e"^ \\9\\l^ 



1 + C 



L2 



Step 5. Proof that 1 + C is surjective on iu„if (M'^). Let 5 G L'^^^f (R'') and consider 
= glp^ for L e 2N+1. As l+C is invertible on L'^{R'^), there exists /l e i^^j^d) 
such that 

(1 + C)fL = 9L (32) 

and 



l/ilL^ ., = sup 



r+i 



1 + . f:^^ 



r+fc5 



< sup ^ 

,2 /cez^'nri. 



1 



1 + /: 



L2 



,-C'(log|j-fc|)2 



Using (Uni), we obtain 
Whhl^^^ < supC ^ e 

for a constant C independent of L. The space L^^^^ (W^) is known to be the dual 

^loc 

(W^), J2kei'' ll/llL2(r+fe) < °°}> which is a separable Banach 
space. Therefore, since the sequence (/l)l>i is bounded in L'^^^f (R'^), there exists 
a subsequence of (/l)l>i (denoted the same for simplicity) and / G ^unif i^"^) such 
that fL / in Ll^^,, (R^) and 



/II,. <liminf||/i||,. _^ <C||.g| 



■^unif /j;_j.oo 



(33) 
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We now want to pass to the limit in the sense of distributions in ([32|). Since C^{M.'^) 
is dense in £^{L'^), the sequence (/l) converges to / in I?'(M'^). Next, we need to 
show that for any ip E 2?(M''), 



(£(/l-/))^-^0. (34) 

We denote by p^^L the density associated with the operator (z — Hq) ^ * {f — 
fL)iz-Hoy\ Then 

jm.<^ ■^I'T^ Jc Jr'' 

and, as ip has compact support, we have by and (pS)) 



where the constant C > is independent of L and z € C. By the dominated 
convergence theorem, it is therefore sufficient, for proving to show that for 

any z G C 



Pz,L'f — > 0. 

! L-4-00 



(35) 



For i? > 1, we define Pz,L,o\it,R and Pz,L,in,R to be the densities associated with the 
operators 



1 



lKd\B(0,fl)^m * (/ - /l) 



1 



and 



^ -ls(o,i?)^m*(/-/L) ^ 



respectively. Therefore p^^^ = Pz,L,out,R + Pz.LAn,R- Let e > 0. In the following, we 
will choose R large enough such that J Pz.L.out.Rf is small for any L. Then, using 
the weak-* convergence of to / we show that J Pz,L;n\.Rf is small for L large 
enough. Reasoning similarly than in the proof of (j22p . we find 



Pz,L,oMt,RV 



<C\\f-fL\\L^ 



(36) 



1 



WL<^\B{0,R.)- 



Now, we need the following lemma. 



WL'^\B{Q,R) 



Z-Hn 



Lemma 3.6. Let W G Ll^^ (M"*) and H = -A + W. There exists C > and 
C > 0, depending only on \\W\\j^2 , such that for any x G L^{R'^) andrj G L°°{W^) 
satisfying R = d {supp{x), supp{r])) > I, and any z G C \ <j{H), we have 

X{z^ H)-' V < Cci(z)e-^'^^(^)« Wvh^ Wxh^ , 

where ci(z) ~ d(z,a{H)y^ , C2{z) — d(z, cr(iJ))/(|z| + 1). In particular, if A is a 
compact set of C \ cr{H), then 



X{z- H) ^77 



where C and C do not depend on z but depend, in general, on A. 
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Proof of Lemma \3.6[ We have 



x{z~H) ^f] 



\x{x)G^{x,y)T]{y)f dxdy, 



where G^{x,y) in the kernel of (z - H) \ By [29, Theorem B.7.2] and [H Corol- 
lary 1] we have for — > 1 

where C > and C" > depend only on ||M^||^2 . Therefore 



x{z^H)-'r, ^ <Cci{zfMl^\\x\\l. sup 

'^2 £c6supp(x) 



□ 



We now go back to (|36p . Using that ||/ — fL\\]^2 < C \\g\\i^2 and Lemma 
we have for R large enough 



Pz,L,ont,RV 



<C^ll5lL^„^^ ll^ll^ie-^'^. (37) 
We can thus choose R such that ([57]) is smaller than e/2. Besides, we have 

/ Pz,L,in.Rf = Tr ( lB(0,_R)>"m * (/ - /i) -TT''' "IT ) 

= / ^BiO,R)Ym * {fh - f)p, 
JS.'i- 

where p is the density associated with the trace class operator (z — Hq)~^ v{z — Hq)~^ . 
For i?' > 0, we have 



'^B(0.,R)Yyn * (/l - f)p 



< 



Ym{x - y) (/ - /l) (y) dyp{x) dx 

(o,fl) 

i^mCa:^ -y){f - fh) (y) dyp{x) dx 



B{0,R.) JB{0,R') 

Ym{--y){f-fL){y)dy 

'\B{0,R') 



(38) 



L°°(B(0,fi)) 



As Ym is exponentially decaying, we can choose R' such that the second term of the 
RHS of (PSI) is smaller that e/4. As to the first term, by the weak-* convergence of 
/z, to / in Ll^.f (R''), we have that 



/il(x) 



Ym{x~y) {f-fL){y)dy 



B{0,R') 



0, 



for any x £ -6(0, R). Besides, we have for a.e. x £ B{0, R) 

\hLix)\ < \\hL\\^^ <C\\f- fLhi^^^ < C \\g\\^.^^^ 

(see psp). By the dominated convergence theorem, it follows that one can choose 
L large enough such that the first term of the RHS of (|38)) is smaller that e/4. 
This concludes the proof of ([35]). thus the proof of (p4| . We are now able to 
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pass to the limit in ((32|) , which concludes the proof of the surjectivity of 1 + £ on 
L^jjjf (R"^). In view of (|33|) . we have shown that there exists C > such that for 
any g £ L^^^.j (M'^), there exists / £ Ll^^^ (M'*) such that 



(l + £)/ = .g and ||/||^._^^^ ^^HffHiL. ' 

Step 6. Proof that 1 + C is injective on Ll^-^^ (R'*). Let / S Ll^^f (If 
(1 + C)f = 0. For B > 1, let xb. = X{o},R in Lemma [3.41 Then, 



(39) 

^) be such that 



and thus 



XRf + XnAf) = 0, 



(1 + £) {xrD = CxrJ - Xi7/:(/) = [C, xr] /■ 



As g :— [C,xr] f G then the solution (p — xrJ of (1 + £)ip = g is unique 

and satisfies 1 1 (^11 ^2 < (^11.911^2 by Therefore 

WxRfhi^^^. <c\\[c,XR]f\\Li^^^ ■ 

Using Lemma 13.41 we have 

WxRfh^ ., <c\\[c,XR]fh^ . 



c 



(40) 



As ||X-r/IIl2 is a non-decreasing function of R converging to ||/||^2 when R — > 
+00 and the RHS of (gO]) goes to when R +00, then ||/||^2 = and / = 0; 
which proves that 1 + £ is injective. The boundedness of 1/(1 + £) then follows 
from . This concludes the proof of Theorem 13.11 □ 



In order to complete the proof of Theorem 13. 1[ we need to prove Lemmas 
andl3.5l 



Proof of Lemma \3.4\ For simplicity, we use the shorthand notation XR = Xi,Ry 
VR = Vi,R and Br = Bi,r. 
Step 1. Proof 0/ (HH). We have 

mf *Y^-Ym* ivRf) = VR (-A + f - (-A + r,Rf 



We now use that [B, {z - Ay^] = {z - Ay^ [B,A] {z - A)''^ and the fact that 
[riR, A] = - {Ar]R + 2Vr]R • V). We thus obtain 

mf * - y„ * {r^Rf) = (-A + to2)-i [^^^ A] (-A + 112^' f 

- - (-A + {{At^r) + 2 (Vr/n) • V) (-A + m^)"! /. (41) 

As VrjR = — Vxfl, and ArjR = —Axr are supported in B2R. \ Br, then, by 



C 

c 



B2r\Br 



lB2«\SHV(-A + m2) V 



(L2(Rd))<* 



(42) 
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To bound the first term of the RHS of ([42]) , we write 

(-A + m^) ^f{x)^i Ym{x-y)f{y)t(^gd\B3R)uBa/.iy)dy 

+ / Y^ix-y)f{y)lB,,\BRjy)dy (43) 



Thanks to the exponential decay of Ym and the fact that for any x G B2R \ Br and 

yeiR''\ Bsr) U Br/2, \x~y\> R/2, we get 

'^B2r\Br (-A + m^) /l(Md\S3„)uBj,/ 



L2 



Yyf * (/1(r<i\b3_„)uBh,2) 



/I 



ControlUng in the same way the second term of the RHS of (j43|) . we deduce 



L2 



< Ce 



l(R''\B3fi)USB,2/ 



C 



We proceed similarly for the second term of the RHS of (j42l) using that = Vym, 
the inverse Fourier transform of i|5'''~^| j^q^^i^j is exponentially decaying and 

satisfies \\Wm * 511^2 < for any g £ . We get 



C R 



H-^^ R 



(44) 



We turn now to estimating ||[?7fl,£] /||^2. We know that [iirjC] f is the density 
associated with the operator 



2i7r z — H{ 
1 

~ 2^77 



(i^ * (j?/?./) - VRYrn * /) ^ dz 



Z-Hn 



1 



Ym*f 



1 



(45) 



We denote by ri and r2 the densities associated with the first and second terms of 
the RHS of (gSl) respectively. For any W E L'^{W'-), we have 









/ riW 




J-/Trf 




2nr Jc \ 



^ ■ (Y„, * imf) - VRYrn * /) -^-ET^ 1 dz 

Z ~ Hn 



(46) 



< C \\Ym * {mf) - VRYm * f\\L2 \\W\\^2 , 
where we have used and (|20p . Therefore, in view of (|44l) . 

Iki||i2 < c||y,„ * (?7k/) - fyK>"m * /IIl2 

c 

1(M''\b3h)uBh/2/ + ;^ 

It remains to estimate r2. For any A e 62 (^^(M'')) and W G L^(]R''), the density 
p associated with the operator (—A + A (—A + satisfies 



C 

< — e -1 
- R 



(47) 



pW 



< 



W\ 



1^1 



e2 



< 



C\\W\\^2\\A\ 



62 
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Therefore 



IIpIIl. <cp|| 



(48) 



Applying (gSD for ^ = (-A + if'^ [r^, {z - H^)-^] Y^*f{z-Ho)-^ (-A + 



we obtain 



z - Hv, 



dz, (49) 



where we have used that Ci(l — A) < |z — i7o| < (^2(1 — A), whose proof is 
similar to the the one of Lemma 13.21 As the commutator [77^, A] has its support 
in B2R \ Br, we consider separately /ls3j,\-Bn/2 ^^'^ f'^iM''\B3B.)vBR/2- Using the 
same techniques as above, we obtain 



(-A + 1)-^ [rjn, A] j^Y,^ * (1s3«\b«/./) (-A + 1)"^ 



< C 



(-A + 1) M^H,A] 



1 



<C{\\Vm\\L- + \\^VR\\L-) 

c 



z-Ho 



62 



L2 



< 



R 



(50) 



Far from the support of [rjR, A], we have 
1 



(-A + 1) M?7fl,A] 



z-Ho 



/) (-A + 1)-^ 



62 



< C 



(-A + 1) MA?7fl,-2V-Vr/H)|| ^ 



B2r\Br 



z-H, 



■Ir+fc 



©2 



lLr+fcy,„ * (l(R^\B3«)ui3«/2/j (-A + 1) ^ 
In dimension d < 3, i/i(M'') ^ L^{W^). Therefore 



(51) 



^\Bsn)UBr, 



./2/ (-A + 1)-^ 



< c 

< c 



tv+kYm * [^(K-i\B3R)yjBn/2 
Ir+fcFm * (l(R''\S3n)US„/2/) 



/)(-A + l) 

< c 



64 

Ym. * \t(m,d 



\B3R)LiBn/2 



(52) 



Using the exponential decay of Ym, we obtain 



Ym * {innd 



AS3r)UBh/2 



\-B3h)U-B^/2 



(53) 



In particular, for £ Z"^ n (i?5_R/2 \ 53_r/4) (the pink part in Figure [3] below) , the 
distance between k and (R'' \ B^r) U i?i^/2 (the blue part in Figure [3]) is greater 
than or equal to i?/4 and 
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Figure 3. Schematic representation of R+ used in the proof of Lemma 



Yrn * ( l(R''\B3J?.)UB„/2/ 



< Ce-'Tre-xd(fc.(R''\-B3R)uBfl/2) 



^{R''\B3R)UBa/2f 

(54) 

Besides, using Lemma 13.61 with i] = £md x = Ir+fc, we obtain 



B2r\Br 



< (^g-C'd(fe,S2J?\-BR) 



In particular for k £ 'Z'^ \ \ ^3^/4) : we have d{k,B2R \ Br) > (see 

Figure [3]) and 



(55) 



Combining dST]), (Ell), (US)), (EH) and ([251), we obtain 



(-A + 1) ^ [rjR, A] ^^i^m * (l(R^\S3«)uB«/,/) (-A + 1)" 



C 

c 

^r' 



'C'R 



'C'\k\ 



l(R<i\B3R)UBH/2/ 



H- 



'C'R 



(56) 



This completes our estimate on r2- Indeed, in view of (|49p . ([SO)) and ([55)) . we 
deduce that 



C ' 

||''2||l2 < — l(E'i\B3H)USi,y2/ 



i? 



1b3h\Bh/2/ 



H- 



which concludes the proof of (P7|) . 

Step 2. Proof of ((28)) . The proof of ((28)) for functions in i^^jf is similar to the 
one of ((27)) for functions. We sketch here the main steps of the proof, and only 
highlighting the differences. Let / € i^nif (K^). Using ((4T)) . we have 

TlRYm * f -Ym* {vrD = E * (Ir+fc/) - i^m * (ryiilr+fc/) 

= ^ (-A + m^)-! ((A77H) + 2(V77fl,) • V) (-A + m^yHr+kf 

kei.'' 

= * (ArjRYm * / + 2Vryfl • VF^ * /) . 

Therefore 

WmYm *f-Y* {VRf)\\Hl^^, < C II (-A + m^) (ry^^r™ ILj^, 

< C II Ary^y™ * / + 2VriR ■ VF™ * /||i2_^^^, 

c c 

- \\^B2R\BRYm * /||i2_^.^ + ||ls2n\BHVy,„ * fW^^^^^^, 



(57) 



22 



SALMA LAHBABI 



To bound the first term of the RHS of ([57|) . we use the exponential decay of Ym, 
the fact that Ym £ f-^iL^) and Lemma [331 We get 



\^B2R\BuYm * f\\^2 . < e" 



< C 



(/I 



l/ll 



h)U-Bh/2 
1b3r\Br/2/ 



As is also exponentially decaying and is in £^{L^), we proceed similarly for 

the second term of the RHS of ((57)) . Finally we obtain the stated inequality 



(58) 



WmYm *f-Ym* < ^ (e-''''' + ||1b3«\s«/,/||^, ^ ) 

We turn to estimating ||[77/;, £] /||^2 . By we have that 
[rjn, C] f = ri + r2 = ri + r2i + ^22 

where ri and r2 are the densities associated with the first and seconds term of (|45p 
respectively, which are now locally trace class operators, and r2i is the density 
associated with the operator 



2l7T Jc 



1 



z-Ho 



1 



z-Ho 



-dz. 



Ym * \^B3a\Bj,/2f) 

By 1231) and using that, in dimension d < 3, iJ^^jf (M'*) L^{W^), we find 

WiWl-^^ < c \\Ym * (W) - mYm * /ILoo < c \\Ym * (mf) - mYm * JWhi^^^ 



-C'R I 



/II 



^B3R\Bfi/2 f 



where we have used (j58p in the last step. Similarly for r2i, since ||A?7fl||^o 
II^'J-rIIl^ — C'/Ri have 



lk2l|lL2_^^j 



< 



[VR, A] —Ym * (1b3«\B«/2/) 



< 



{(Am) + 2{Vm) ■ V) 



1 



Z-Ho 



B 

Ym. * 



(li33B\Bfl/2/) 



< 



c_ 

R 



'^b^rXBr/^I 



(59) 



As to r22; it is actually in L (M ) and 



|?-22||i2 < ||r22||£2 < 



1 



(-A + 1) [r^R, A] —rrYm * (1(r.\B3«)ub«/,/) 



< 



/II 



z-Ho 

'^Bs.rXBr/^I 



S2 
(60) 



The proof of (j60| is exactly the same than the proof of ([56| , except that in (|53p , we 
use the inequality \\Ym * /H^oo < C 11/11^2 instead of the inequality \\Ym * JWh^ — 
C \\f\\ff-i ■ This concludes the proof of the lemma. □ 



We pass now to the proof of Lemma 13.51 
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Proof of Lemma lKR We denote by ?/„ = Xa" and 6„ ~ Ca^"^e^^ " for rt G N and 
a > ao ~ max{a, 2}. By the assumption ([29]), /a < Xcn-i = and we have 

C , C 

Vn < Ona;o + —Xar^/a < hnXQ + — yn-L 

a a 
Besides, there exists a continuous function C(a) such that < C(a)z„, where z„ = 
(7"/q,"("+i)/22;(j is a sequence defined by the induction relation z„ — C/q!"z„_i. 
Going back to xr, we deduce that for any n G N \ and R = a"', we have 

log(i?.)'' I ^11 los(H) ^' l°g(«)^ , , 

a;fl < C(a)e~^ -n^i(^+<^ t^^xq < C(a)e"^ "f^^iM-xo- (61) 

As (pT|l holds true for any a G [ag, ag] , we deduce that there exists C > inde- 
pendent of a, but depending in general on a, such that for any i? > 2, 

which concludes the proof of the lemma. 

□ 

4. Proof of Theorem 12.11 (Existence of ground states) 

Let us now establish the existence of a ground state for the perturbed crystal in 
the rHF framework. The proof of Theorem 12.11 is a consequence of our results on 
the operator £ stated in the last section, and of the properties of the higher-order 
term in the expansion of for a charge distribution / G iynif (M'*). 

To solve the self-consistent equation (jTTI) , we first formulate the system in terms 
of the response electronic density p = py — p^^ as follow 

PQ 

Q = 1ho+v„<o — 1ho<o (62) 

-AK + mV, ^ IS-^-^ {p-v). 

Indeed, if p is solution of (|62|) . then 7 = 1 (i/o + Ym * {p ^ Pjo ~ i^) < 0) solves ([Tl]l . 
For a charge density / G L'^^-^f (K'^), we expand 

Qf = liHo+Yr^*f<0)~^l{Ho<0) 

as powers of / when / is small. For this purpose, we assume that 

d{C,a{Ho))>g, 

where g = d(0, a^Ho)) and C is now a smooth curve in the complex plane enclosing 
the whole spectrum of Hq below and crossing the real line at and at some point 
c < inia{Ho) - g (see Figure[2]). Let us recall that for V G L°°{R'^), a {Hq + V) C 
a (Ho) + [- ll^llioo , \\V\\^^]. Therefore if < g, then Ho + V has a gap 

around and a (H) C [inf a {Hq) — g, +00). For such a V, we have using Cauchy's 
residue formula, 

Q = 1 (i/o + F < 0) - 1 (i/o < 0) = / -J -dz - / -^-TT^^^- 

2i7r Jc z- Ho-V 2nT Jc z- Hq 

By the resolvent formula, we obtain 

' ' ' V^dz+^U^vY d.. 



2iTT Jc z ~ Hq z — Hq 2iTT Jc \z — Hq J z — Hq — V 
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Therefore for / e Ll^^ (R"*) such that ||/ * Y^nho. < g, 

Qf = Qij + Q2j, (63) 
where Qi.f has been defined and studied in Section [3] and Q2,f is defined by 



2i7r 



z- H, 



■Y,n * f 



z- Ho- Y„i * f 



dz. 



We give some properties of the second order term Q2J in Lemma l¥!T] below. Using 
the decomposition (|63p . equation (p^ becomes 



P = PQl.p-„ + PQ, 

FoUowing ideas of [15], we recast as 



= -i^{p -^)+ Pq^ 



C 



P = 



1 + C 



1 + /:^Q2(p-f)- 



(64) 



(65) 



In Proposition 14.21 below, we show that for ly small enough, the operator : 
p I— > £ (1 + £) ^ + (1 + C)~^ PQ2(p-iy) ^"inrit^ ^ fixed point, which is controlled 
in the i^jjjf norm by the nuclear perturbation v. This will conclude the proof of 
Theorem [O 



Lemma 4.1 (Properties of the second order term). There exists > and C > 
such that for any f £ L^^if (M'') satisfying \\f\\]^2 < Sc, the operator Q2J is trace 
class, the density is in i^^jf (W^) and 



<c\\f\\l. 



Proof. Since \\Ym * /H^oo < Co |l/||i2 ^ (see ([25])), we can choose 5c = g/2Co, 
where, we recall that g — d(0, cr(i?o))- In this case, {z ~ H ~ y,„ * /)^^(— A + 1) 
and its inverse are uniformly bounded w.r.t z e C (see Lemma l3.2p . Using the 



exact same procedure as in the proof of (|231) . we obtain that Q2.f is trace class, 
Pq.j e Ll^^i (J^") and 





< C 







Yrn * f^-FT 

c z- Ha 



Y„i * fdz 



<C\\Y 



\\Y,r.*f\\U<C\\f\\l.^^^ 



which concludes the proof of the lemma. 

Proposition 4.2. There exists Q!c,e > such that if ||j^||^2 < ac, then 



□ 



P ^ 



Bl^- (e) 



l+£' 



is well-defined and contracting onBj^-i ^ (e) = |/ G -^unif (1^'^), 11/11^2 — ^^'^■'^ 
it admits a unique fixed point p in the ball Bj^2 ^ (e). Moreover p satisfies 



(66) 



for a constant C independent of v. 
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Proof. We want to use Lemma |4TT] to show that Q is well-defined on a small ball of 
L^jjjf (R"^). Here, the charge distribution \s f = p — v. We thus need to choose Uc 
and e such that \\p — r2 < ||p|| r2 + r2 < e + ckc < (^c, where 5c is given 

^ unif unif unif 

by Lemma |47T] Let yl > 0, < e < 5c/ {1 + A) and ac = Ac. Let and p such that 
II'^IIls < ac and ||/9||^2 < e. By Lemma [47l1 and the fact that £ and 1/(1 + C) 

unif unif 

are bounded on L^^^-^f (U.'^) (see Theorem 13. ip . we have 



\5M\ 



< 



c 


Ml- ,. + 

6(iL.t ) 


1 








1 + C 


1 + C 


BiL^,, ) 







<C'ilkllL;„j +C2||p-t^|li2_^^^ < {CiA + C2il + A)h)e. (67) 
We choose A < 1/Ci such that for e < (1 - ACi)/(C2(l + ^)^), we have 



< e. 



To show that Qi, is contracting on i?^2 ^ (e) for e small enough, we use the explicit 
expression of Q2.p-u- Let p, p' G Bip. ^ (e) and denote by _ff = + y„j *(/?—;/) 
and H' = Hg + Y„i * (/?' — v). The function (1 + £) {Gu{p) — Gv{p')) is the density 
associated with the operator 



2iTT Jc\Z~ Hq 



Ym *{p~- v) 



z- H 



1 



z-Ho 



Y,n * {p' - v) 



z- H' 



■ dz. 



A straightforward calculation shows that this operator can be written as 

2 



1 



1 



Ym *{p-v) 



Yjn*{p~ P') 



z-H 



z-H' 



1 , , 1 . A 1 

Yrn*[p- V) ET-^™ * [P - P) 



+ 



z^H„ 
1 



Z-Hn 



Ym*{p- p')- 



-Ym * [p' - v) 



Z-H' 
1 



z-H' 



dz. 



(68) 



2; - iJo ' Z- Hq 

Using the same techniques as before, we deduce that 

WQuip) - Qu{p')\\l-- < C3 (||Hli2 + \\p'h. + ||Hli2 )\\p- p'h. 

unif \ unit unit unit / unit 

<C3i2 + A)e\\p-p'\\^.^^^ . 

Taking, in addition, e < 1 / {C3{2 + A)) , we have that Q,y is contracting on Bj^2 ^ (e). 
Let p be the unique fixed point of Gi, in _B^2 ^ (e) . It remains to prove (|66|1 . By (|67p , 
we have 



IIpII 



\\GAp)\\l- , <C^ilkllL^ , +C2(l + A)e(||p||^2 +||Hli2 ) 

unit unit \ unit unit / 



Therefore (1 - 6*2(1 + A)e) \\p\\^2 < (Ci + C2(l + A)e) \\i^\\^2 . Using that e < 
(1 - ACi) I (C2(l + Af), we have 1 - C2(l + yl)e > and we"deduce that 



^ Cl+C2(l+A)£ 



1-C2(l + A)e 
which concludes the proof of the proposition. 



< 



A 



□ 
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5. Proofs of Theorem 12.31 and Proposition 12.51 (Decay estimates) 



We present in this section the proofs of Theorem 12.31 and Proposition 12.51 They 
consist in decay estimates of the mean-field potential Vi, and the mean-field density 
p,y. These estimates are used later on in the proofs of Theorems 12.21 and 12 . 71 



5.1. Proof of Theorem [im 

Proof of Theorem ] 2. 31 Assume that 



1^2 < ttc, where ac is given in Theo- 
rem [^H] We use the notation p to denote the mean-field density p^, — p^^^-^g, the 
solution of and denote hyV — V,y — Y„i*{p—i'). Recall the decomposition ([M)) 
of p in a linear term and a higher order term 

p=-Cip~iy)+PQ^^_^. 

Using localizing functions, we will show that p decays far from the support of v. 
To do so, let us introduce the set / = {fc € Z'^, supp(i^) n -8(0, 1) -I- fc 7^ 0} and for 
i? > 1, the set Bfi = Bj^ji and the the function XR = Xi,R defined in Lemma 13.41 
We denote by ry/f = 1 — XR- We thus have 

'HRP = -^VR^iP -1^)+ VrPq^^_^ = -^VRiP + [C, Vr] {P~'')+ VRPQ^^_y 
As for i? > 1, rjiiiy — 0, it follows 



1 



VRP = 



(1 + ^) 



[C, r]n] [p - v) 



1 



(69) 



We will successively bound each term of the RHS of ([B^ . For the first term, we 
have by Lemma [3.41 for i? > 2, 



1 



[C, tib] {p - v) 



C 

< — 

- R 

C 

< — 

- R 



e-^'^\\p-v\\^. 



lB3«\i3R/2 ip-^) 



^B3r\Br/2P 



(70) 



where we have used that 1^3/1X3^/2'^ ^ ^ — 2, that p is controlled by v in the 

i^jjif norm and that 1/(1 +£) is bounded on L^^^^f iW^). As to the second term of 
the RHS of since l^d^s^Tyi} = rjn, we have 



1 



2iTT 



Ho' 
1 



1 



z-Ho 



1 

2i7r jQ 



VR, 



z - 
V- 



Ho 

VR, 
1 



VrV- 



1 



z - 



-V- 



1 



H 
V- 



-dz 



1 







H 



z- H 
-dz. 



-dz 



(71) 



where H = Hq + V and C is as in the previous section. We recall that by the 
assumption ||i^||^2 < ac, the operator H has a gap around 0, thus the operator 
{z — H)~^{—A + l) and its inverse are uniformly bounded on C and all the estimates 
obtained in the previous sections hold when we replace Hq by H. We denote by 
rs, r4 and rs the densities associated with the three operators of the RHS of ([7T|l 
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respectively. Using an inequality similar to (|23p . involving H instead of Hq in the 
resolvent in the right, we have 



1 



V 

" z-Ho 



VriR 



dz<c\\viR.\Bnh^ WVmh^ 



By ([28| in Lemma [3.41 and using that \\Y„i * /||^2 — 11/11^2 , we have that for 
R>2 



B3r\Bii/2P 



(72) 



Therefore 

1^311^2 <C||1 



-C'R I 



C_ 

R 



To bound r4 and r^, we recall that we have shown in the proof of (|28p (see ([5^ 
and (HOI)) that for any / S L^^jj (R'') 



(-A + 1) M'7i?,A] 



z- H< 



■Yrn * (1(R<'\B3r)UB2h/) 



■Y,n * (1 



B3r\B2R 



/) 



< 



62 

-C'R 



11/11 



^3R\B[i/2.f 



Therefore, using again the equality \r|B.^ {z — Hq) ^] = -~{z — Hq) ^ [r//j, A] (z 
Ho)~^, and an inequality similar to ([23| . we obtain that for any R> 2, 



<^l|l 



Z-Hn 



-V 



dz 



B3r\Bji/2P 



(73) 



The last term of the RHS of (1711) can be written Qia + Qout , where 



1 

2j7r ./c 



1 



z-Ho 

In the same way we obtained ([75)) . we get 



* (li33fl\S2„ (P - J^)) 



B-3r\Br,2P 



To estimate pq^^^, we recall that by (|48)) . we have that for any A E 62(^^(1^'*)) 



Therefore 



iPQautllLl^.^ 



< c 



(-A+1)" 2 A(-A+l) 2 



(-A + l)-M^fl,A] 
1 



L2 



1 



y(i-A)- 



dz 



<C\\A\\ 



-Ym. * (IcKd 



\B3h)UB2 



B3r\Bj^ /2P 
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Now that we have found estimates on rs , r4 and = pq^^ + pg^^j , we use that 
to estimate VRPq^ follow 



B3r\Bii/2P 



(74) 



Using once more that 1/(1 + C) is bounded on L^^jf (R^), we deduce in view 
of (ESI), (EDI) and (El 



B3r\Br/2P 



L2 



We choose < min {l/(2Co), ad, and assume that ||i^||^2 < a^. It follows 



C 



1/ 



We have a similar inequality for V. Indeed, by ([7^ . we have 



unit un 



c 










C 







c 



Using Lemma 13.51 with xjj to ||l]gd^^^p||^2 , we obtain 



(75) 



(76) 



□ 



Inserting (1751) in ([75)) . we get 

Finally, noticing that lR<i\Ci?.(i^) — '7^/2, we conclude the proof of (1131) . 

We now turn to the 
5.2. Proof of Proposition [2751 

Proof of Provosition \2.5[ Assume that || 1^11^2 < ttc, where etc is given in Theo- 
rem As pi, and pi,j^ are fixed points of the functionals Qi, and Q^,^ respectively, 
then 



C 



Pu - Pui^ = 



1 + C 



(i^ - vl) + 



1 



1 + £ Q2(p,.-l^)-Q2(p,.i-l^L)' 



For i? > 1, let xr — X{a},R Bji — B{g} as defined in Lemma 13.41 Since 
1b„ < Xfl, then 



< 
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Besides, there exists C > such that for any / £ L^nii i^''') and any R> 1, 



29 



1 



XR, 



1 + C 



I 



4 



11/11 



1b3 



1 



l + C 



(77) 



Indeed, using that l/(l-|-£) is bounded on i^^jf (W^) and estimate (PSI) in Lemma 
(notice that Lxr — XbJ^ — VbJ^ ~ ^Vb)i '^^ obtain 



1 + C 



f 



< c 



1 



< — e 
- R 



C / 
< — 

- R 



-C'R 



Li, 







' f 

1 + C-' 


+ 

-'"unit 



^'^11/11 



Using dZZl) for / = £ (i^ - vl) + 



we have 



1 + /: 



XRCiv-VL) 



1 



1 + /: 



, +§I!1b3h(p.-p.JIIl^ , ■ (78) 

We first bound the first term of the RHS of Using (gH) in Lemma [23 and 

that for i? < L/4 it holds that Xfl(i' — (^l) = 0, we have for R < L/A 



l + C 



XrC{v~vl) 



<C\\xRC{u-iyL)\\L- 



c 



C 



C 



-C'R 



c 



(79) 



We turn to tlie second term of the RHS of ([78]) . Using ([28)) . a decomposition similar 
to ([55)1 and reasoning as in the proof of ([73]), we find for R< L/A 



1 



^ ^XRPq^(^P^_^)_q^(^P^^^_^^) 



/ C ' c 

We choose < Uc such that 6*0(0:^ + ct'^^) < 1/2. Thus, if ||i^||^2 < then 



1 + C 



XRPr 



+ §(^""'"llHl.L.+l|l-3./.ll..„,). (80) 
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In this case, combining ([78l) . ([79|) and (|80|) . we obtain for i? < L/4 

Using a recursion argument, we easily see that for any /3 > 1, there exists C > 
such that 

C / C 

c c 

To conclude the proof of the proposition, it remains to prove the bound on the 
potential. Using (|28|) and denoting by Jl = Pv — Pu^ — ^ + ^l^ we have 

IIK - KJIh2_^^^ {B^,^p) - \\XL/Af<Yir, * IlWhI^^^ 



L 



-C'L 



C 



C 



□ 



6. Proof of Theorem 12.21 (Thermodynamic limit) 
Proof of Theorem \2.2[ Assume that ||i^||2^2 < ac, where etc is given by Proposi- 



tion [2T5l By Cauchy's formula, we have 



1 



7u - lui^ = 



1 



1 



-dz, 



2iTT Jc z — Hq — Vy z - Ho — Vu^ 
where the curve C is as in Section S) We write the resolvent difference as 
1 1 1 _ „ 1 



z-Ho-Vu z-Ho-V^^ z-Ha-Vu - z-Ho-V^^ 
where f l = Pv — v — p^^ + v^. For a compact set _B C K'', we have 



Tr 



1b 



1 



z - i7o - K 



* fh 



1 



-1; 



< C 



1/ 



1 



* fh 



For L large enough, we have B C -6(0, L/S) and, by Proposition 



1b 



1 



z-Hq-V, 



ls(0,-L/4)^)n * /i 



< 



Is 



1 



|lB(0,L/4)i^m * fh 



c 



<^\\v\ 



Besides, as d(i3,R'^ \ > L/8, we have using Lemma 



z-H„~V, 



■lR'i\B(0,L/4)^m * /l 



< Ce 



-C'L 



lK<i\B(0,L/4)^m, * /i | 
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As C is a compact set and all the estimates are uniform on C, we conclude that 



0. 



□ 



7. Proof of Theorem 12.71 (Expansion of the density of states) 



The proof of Theorem 12.71 follows essentially the proof of [19, Theorem 1.1]. The 
main difference is the proof of Proposition l7.1l below. which deals with self-consistent 
potentials, while [191 Proposition 2.1] deals with linear potentials. Treating non- 
linear potentials is done at the price of assuming that the defect x is small in 
the iynif -norm, so that the potential decays fast enough. For the sake of self- 
containment, we mention here the main steps of the proof; more details can be 
found in [21]. 

Proof of Theorem \2.7\ . Following [T9] , we first express the density of states of the 
random operator Hp{uj) in terms of the resolvent {z — Hp)^^ for z G C We 
next find an asymptotic expansion of Tr ((z — Hp)^^^ using a thermodynamic limit 
procedure. 

We recall the Helffer-Sjostrand formula ^18. .llj. For a self-adjoint operator A 
and if E 5(M), we have 



ip{A) = -- I ^{z)^—rdxdy, 
IT OZ z — A 



where : C — )• C is an appropriate complex extension of ip such that 

(i) (^e e C, |Im(z)| < 1}), 

(ii) for any n G N and a, /3 > 0, one has 



sup Na,p \Xi-^ \ \y\ 

|y|<i 



dz 



[x + iy) ] ] < C„^a,p sup Na',l3'{^), 



3'<ii + /J + 2 



(81) 



where Na,0{(p) = sup^^j 
Hence, for ip e 5(M), 

{np-no,(p) = / (p{x)np{dx) 



^ dxf 



TT \Jc dz 

Besides, denoting by = Vj^^ , we have 

1 1 



ip{x)no{dx) = Tr {p{Hp) 
dx dy 



z - H„ 



z - H.„ Z - Hn 



Z~Hn 



Therefore, using the Kato-Seiler-Simon inequality (|19p and Lemma 



Tr 



Z — Hr, Z — Hn 



< 



lr(-A + l) 
1 



(-A + 1) 



Z-Hn 



we obtain 

^llL°°(nxR<i) 



(-A + i; 



(-A + l)"'lr 



< C 



fl + \z\_ 
{\lm{z)\ 



PllL°=(n> 
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By Fubini's theorem, we get 



{up -no,ip) 



dy. 



(82) 



In the followhig, we find the asymptotic expansion of Tr {{z — Hp)~^ — (z — iJg)^^) 
asp — > for z G {C \ R, |Im(z)| < 1}. To use a thermodynamic hmit procedure, we 
consider, for each reahzation cj S 17 and each box size L E 2N + 1, the system with 
defects only in the box F^, that is, we consider the defect distribution VKLiu:)i^)-< 
with Kl{uj) = {keZ'^n Tl, qk{uj) = l}. For K C Z'', we recall the notation 
'^K = Y.k£K X(- -k),VK ^ VvK ^Ym* {Pi^K - ^k) and Hk ^ Hq + Vk ■ By the 
proof of Theorem 12.21 we have, almost surely. 



Tr Ir 



Besides, from (|T2j) and (|25]), it fohows 
1 1 



Tr Ir 



z-Hp{uj) z-Hk^(u) 



li 



< C 



VIH {z) 



\\Vp{lo,-)-VkA^,-)\\ 



The dominated converge theorem thus gives 



and 



Tr 



z — Hp z — Ho 



^ = lim E I Tr ( Ir ( — 

L^oo \ V \z-Hkl z - Hq 



(83) 



Let L e 2N + 1 and N = L'^. As the random variable Tr(lr((z - Hk^)^^ - 
(z — i7o)^^)lr) depends only on the N independent Bernoulli random variables 
{qk)kei.'' nrr.^ we have 



E Tr Ir 



z - Hkl z- Ho 



xcz-inrL 

N 



1 



\K\=r, 



z — Hk z — Ho 
1 



z - Hk z- Ho 
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Expanding the term (1 —p) " as powers of p and rearranging the sums, we obtain 

j=0 Ti=0 ^ Kc-L'^nv^ \ \ I' 

\K\=r, 

,7 

= Y.''3^LP' +RjAz,p). (84) 
where we have denoted the j'^ order term by 

-M-t E Mr"Tr(lr(^- 

\K\=-n. 



and the remainder of the series by 

N 



i?,.,(z,p) = 5:p"(l-pr-" ^ Tr(lr(— ^--^)lr)-E«.. 



LfC| = n 



The result wiU now fohow from the next proposition, whose proof is postponed 
until the end of the proof of the theorem. 

Proposition 7.1 (Estimates on a^^^ and Rj^l)- There exists a,, > such that 
• for j < 2, there exists C > such that for any x S L'^{M.'^) satisfying supp(x) C F 

and 11x11 L2 < etc o,nd any z G C \ M, 



K'CK 



V|Im(z)| 



i?,7,L(z,p)|<C||x|lr^P'^+^f^i^) . (86) 



(85) 

• for J < 2, there exists C > swc/i that for any x G -/^^(K'') satisfying supp(x) C F 
and ||xIIl2 < z e C \ R, p G [0, 1] and L G 2N + 1 

, (.7+2)(d+l) 

V|Im(z)|, 

We deduce from Proposition 17.11 that for any j < 2, and z G C \ R, ajx{z) 
converges as L — cx) to 

= E E - . (87, 

and that for any J < 2 and p G [0, 1], Rj^l{z,p) converges, up to extraction, as 
L — 7- oo to Rj{z,p), which satisfies 

Il2 



\RAz,p)\<C\\xh.p-'+' 



fl + \z\_ 
V|Im(z)| 
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Passing to the limit as L — >■ cx) for this subsequence in ((84l) and in view of ([83|). we 
obtain 

\z-Hp z-HaJ ^ 
Going back to ([5^ . we thus have 

A simple calculation shows that 



«.W-7 E E E (-l)l"^"'lTr(lrM 



/f|=3, oeff 



1 



7 E E (-1)1"^"'!^ M 



z - Hk' z - Ha 
1 



It 



J i^^!^ V ^ ^ z- Ha, 

Kczt* K'CK ^ ' 

\K\=j, 0€K 

Therefore, by the dominated convergence theorem for series, we obtain 

/ ^iz)a,iz)dxdy=- J2 E (-l)'''^''''Tr(^(ifKO-^(i?o)) 

\K\=j, oeK 

Moreover, using (ISTjl . (|85|) and ([87| . we see that /ij is a distribution of order at 
most J ' + 3 + jd. Finally, (p i— > —-^ /j, ^{z)Rj{z,p) dxdy defines a distribution of 
order at most J + 4+ {J + 2)d and satisfies 



- / ^{z)Rj{z,p)dxdy 
IT Jc oz 



< Cj sup Na,p{^). 



/3< J+4+(J+2)d 

»<(J + 3)(£i+l) 



This concludes the proof of Theorem 12.71 



□ 



To complete the proof of Theorem 12. 7[ we need to prove Proposition 17.11 We 
first state and prove Lemma l7 . 2 1 which will be useful in the proof of Proposition [73] 



Lemma 7.2. Let H = -l\ + W, with W G L^^jj {W^). Then, for any (3 en and 

any Borel set B C , there exists C > and C > such that for any z G C \ M 
and any v,v' G L'j.iW^) satisfying \\v\\^2 , ||i^'|1^2 < ac, R = d(supp(z^), 0) > 1, 
R' = d(supp(i/'), 0) >1,D = d(supp(z^), supp(!^')) > 1, we have 



z-H 



|Im(z)| 



1/1 f 



z-H 



( W - K) 



< 



c / 1 



|Im(z)| V^*^ 



(89) 
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(90) 



and 



z- H 



< 



iim(z)i ^i?'^ 



(91) 



where R = min{i?, i?'}, C2(z) = A{z,<j{H)) j [1 + |z|) and where the constants C 
and C" depend on W only through its L^^^^f -norm. 

Proof. Inequalities - (PT|) follow from Lemmas 13.21 and 13.61 Theorem 12.31 and 
Proposition 12.61 For instance, for ([55)l . we first look at Vi, far from F. Using 
Lemma 13.61 we have 



It 



\B(0,f ) 



< 



|Im(^)| 

C 
|Im(z)| 



Near F, K decays as R gets large by Theorem^ As d{B{0, f ), supp(i^)) > i?/2, 
then, by ([T5)) . we have 



S(0,f) 



< c 



S(0.f ) 



<^g-C'(logfl)^||^|| 



where we have used that in dimension c? < 3, H^j^^f (M'') ^ L°°(R''). We next use 
Lemma l3T2l and the Kato-Seiler-Simon inequality (|19|) to obtain 



It 



-1 



z- H 



B(0,f) 









< 






©2 


-A + 1 


62 



(-A + 1) 



H 



B(0,f ) 



< c 



|Im(z) 



g~C'(logfl)^ n II 



which concludes the proof of ([55)) . The proofs of , ([M)) and (PT|) use the same 
techniques; they can be found in J^I\. □ 

We now prove Proposition 17. II 

Proof of Proposition I7.il Let ac be the minimum of the constants Oc defined in 
Theorems 12.11 and 12.31 and Propositions 12.51 and 12.61 We assume that ||xIIl2 < etc. 
Throughout the proof, /3 will denote an integer greater than d+1 whose value might 
change from one line to another and C > and C" > constants that depend, in 
general, on /3. For z e C \ R, we denote by Ro{z) = {z — Hq)~^ and for any 
K C 'Z'^, we set Rk{z) — {z — Hk)~^ ■ We omit the dependence on z when there is 
no ambiguity. We also omit the HxIIl^ estimates. Let L E 2N+ 1 and denote 

by iV = L"^. 

For j ~ 1 and K — {k}, with k E Z,'^, we have 
|Tr (Ir {R^k} - Ro) tr) \ = |Tr (lri?o^{fc}i?{fc}lr) | < ||lri?o^{fc} ||s, ||^{fc}lr||e. 
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Therefore, using ([88]) in Lemma 17.21 we get 



|Tr {lrR,.,V,,,Rolr)\ < C (l±M)' (e-^Co^l'^D ^ e"^- ^I'^l) . 

Since the series '^keZ'' ■ with A > 0, is equivalent to f^^ e^^^^^dx = l/A'', 

and for z e {z £ C, |Im(z)| < 1}, it holds l/c2(z) < (1 + \z\)/ |Im(z)| and 1 < (1 + 
\z\)/ |Im(z)|, we deduce that the series X^fcezt* 1'^''' {^rR{k}V{k}Ro^r) \ is convergent 
and its sum satisfies 

E (lri?u.}V{,}i?olr)| < C ^i±M.^'^\ 

For j ^ 2 and K = {k,k'}, with k,k' £ Z"^, a straightforward calculation gives 
^ (-l)l^\^'lTr(lr(i?,^,-i?o)lr) 

K'CK 

- Tr {trRo{Vik,k'} - V^k} ~ V^k'})R{k.k'}tr) 
+ Tr {trRoV[k}R{k} {^{k.k'} " ^{fe}) -R{fe,fe'}lr) 
+ Tr {trRoVik'}R{k'} iV{k,k'} - i?{fe,fe'}lr) • (92) 

Using the inequality (^11) , the first term of the RHS of ([M)) can be estimated by 
||lri?o(V{fe,fc'} - T/{fc} - V{fc'})||g^ ||i?{fe,fe'}lr||g^ 



< 



f 1±MV I ^ f -C'(logmin{|fe|,|fc'|})' -C'c2(z)min{|fc|,|fc'|}\ 

\\lmiz)\J \ \k-k'fy ) 



V|Im(z)|; \\k-k't 
As to bound the second term of the RHS of ([M)) . it is bounded by 
||lr^oV{fc}||e^ ||^{fc} (V{fe,fe'} - ^{fe})||e ll^{fe:fe'}lr||e^ 

V|Im(2:)|y V J 

using and by 



< C 

using Therefore 



|Tr (lri?oF(fe}i?(fe} (n^,^'} - ^{n) ^{^.^'}lr)| < C 

We have the same bound for the third term of the RHS of (j92|) . Therefore, the 
series Ekc- EK'cK(-l)''^^''"'Tr (Ir [Rk' - i?o)lr 

IK\=2 



is convergent and its sum 



SHORT-RANGE QUANTUM CRYSTALS WITH DEFECTS 



37 



satisfies 



E 

|K|=2 



^ (-l)I^WlTr (Ir - i?o) Ir; 



K'CK 



< C 



(Mil 

y\lm{z)\ 



3+2d 



We turn to the proof of the estimate on the remainder Let J < 2 and 

p £ [0, 1]. We first write Rj^l{z,p) in the form of the expectancy of a binomial 
variable. Indeed, we have 

N J 



n=0 



i=i 



where fL,K — Tr (Ir — i?o) Ir)- Rearranging all the terms (see [21] for details), 
we obtain 

= p-^+iE (gj^L (Fl + J + 1, z)) , 

where Yl is a random variable of binomial distribution of parameters p and N—J— 1 
and gj,L{-, z) : {J + 1, ■ ■ ■ , N} R is defined by 

/ /V — ./ — I \ ^ 

9J,L{n, z) 



N- J-l 
n- J-l 



TT{lr{RK{z)-Ro{z))lr) 

\K\=r^ 

E E E (-l)^''^"Tr(lr(i?K"W-i?o)lr) 



KCZ'^nTj^ K'CK K"CK' 
|K|=„ \K'\<J 



N- J-l 
n — J 



In order to prove ((86| , it is therefore sufficient to show that there exists C > such 
that for any L G 2N + 1 and J + 1 < n < TV, 

It is sufficient to prove the above inequality for J = 2. Let J + 1 < n < iV 
and consider a configuration K C Z'^ HTl such that \K\ — n. A straightforward 
calculation shows that 

/ /V — ./ — I \ ^ 

9J,Lin,z) 



TV- J- 1 
n- J -1 



where 



and 



^ Tr (lri?o iPi,K - P2,k) RkIt) 

\K\ = n+J + l 



k^K 



k,k' eK 



P2.K = ^w^f^ (^^ - ^{fe}) + E {y{Kk'}R{k,k'} {vk ~ V{k.t 

k£K {k,k'}cK 

- y{k}R{k} {vk - V{k}) - y{k'}R{k'} {vk ~ v^k'j 

Besides 

Pi,K = Y ^r+rPi,K. 



'}) 
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For each r £ Z'', we split 'ir+rPi,K into two r-dependent quantities: a part involving 
the defect in fco — arginikeK \ k — r\ and the rest. We denote by 

Arm =Vk - V{ko} - X! (^{'=.'=0} ~ ^{ko}) 
keK\{ko} 

and 

BkM = X! i^{k:k'} ~ V{k} - ^{fe'}) • 
{k.k'}cK\{k„} 

Then 

Pl,K — ^ lr+r^X,fco — ^ lr+r-B_ff,fco- 

We have thus split gj,L{n, z) into three parts 



9J,L{n,z) = M ^ Tr [ lri?o X! ^t+tAkmRk'^t 



N - J -1 
71 - J - 1 

iV - J - 1 
n- J- 1 



^ Tr ^ lr+rB;f,fe„i?A'lr 

|A'|=„ 

^ Tr iP2MRKlr) (93) 



IK 1=7. 



that we will bound successively. We start by the first term. Let r e Z'* and denote 
by fci = arginffcg/^\{fcjj} d (fc, {r, fco})- We introduce 

iQiK,r) - |r-fco|, £i{K,r)^diK\{ko},{r,ko}) 

and 

£2(if,r) =d(if \{fco,fci},{r,fco,fci}). 

When there is no ambiguity, we omit to note the dependence of these quantities on 
K and r. By Theorem 12.31 we first have 

(,to + J- J 

We now want to control || Ir+r {Vk - V{ko}) by l/(£i + 1)'^. If ^ < i^/Af^ (see 
Figure 3]), then by Proposition 12.61 we have 

l|lr+.(^K-nM)IL==< (95) 
If £o > ^1/4*^, then dH]) gives 

||lr.,(V'K->',M)L.<^<^^ 
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Therefore, by ([M)). ([551 and ([M]) 

\\lr+r{VK-V[,„y)\\^^ = 

< 



Ir+r {VK-V[k„})\\l^ X ||lr+. {VK-V[k,})\\l 
C 



(4 + 1)^ (4 + 1)^ 



(97) 



We proceed similarly for the remaining term of AK,ko- First, as (|97p holds for any 
/? > and any K 9 fco, then we have for any k G K \ {ko} 

C 



|lr+r (^{feo.fe} - ^{feo})| 



< 



Next, if i'o < |fc — fcol , then by Proposition 12.61 we have 



(98) 



|lr+r {V{ko,k} - V{ko})\ 



< 



c 



13' 



Otherwise, by 



|lr+r {V{k„^k} - "^{fco})|L=o < 



C 



< 



C 



1^°- ^ (^0 + 1)^ - (|A:-fco| + l)'^' 
Therefore, reasoning as in (|97|) . we have for /3 large enough 



Ir+r {V{ko,k} - y{ko}) 

keK\{ko} 



< E (llw(nM.fc}-%o})ll! 

keK\{ko} 



< 



< 



(€l + l)^(4 + l)',,^fc„} \ko~k\' 

C 

(€i + 1)4(4 + 1)^' 



(99) 



As and ([7f|) holds for any /3 > 0, then by the definition of Akmo, obtain 

C 



(^+1)^(4 + 1)^- 

To control Ax^kn by 1/^2 1 '^^ rearrange the terms of Akmo follows 

^K,fco =Vk ~ V{koM} - E (^{fc^fco} - ^{fco}) • 

i<"\{fco,fci} 



(100) 
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By Proposition 12. 6| we thus have 



< 



c 



E 



c c 

< 



(101) 



As ()100p and (jlOip hold for any (3, then reasoning as in the proof of (jOT)) we have 

Tr (|lr-Rolr+r^if,fco-R_ff lr|) < lllr-Rolr+rlle^ l|lr+r^A:,fco 11^== ll-^-fs'lrllg^ 



< C- 



-C'c2iz}\r\ 



l + \z\ 



IM^)I (£o + 1)^(^1 + 1)^^^ llmWI 



Therefore J^reZ''^^ (llr^olr+r^K./co^iflrl) is a convergent series. By Fubini's 
Theorem, we thus have 



J2 Tr lri?o J2 ^r+rAK,ko RrIf = E E (lri?olr+r^if,fco^if Ir) ■ 



K| = n 



To perform the sum over the configurations K e {^K <Z 'L'^ C\T i^, \K\—n^, we 
classify these configurations depending on the value of li{r, K). i e {0, 1, 2}: 



J2 Tr I lri?o I J2 ^r+rAKM ] RrIt 



VdL 

E E 

r£l'^ Lo,Lx,L2=ti Kci^'^-nr^, IK\ = 



-C'c2{z)\r\ ( 1 + \A 



■/dL 



< f-, ~C'c2{z)\r\ ( 1 + kl \ NL,n,r{Lo,Li,L2) 

^ \\\Mz)\) (L, + ' 

where NL^n,r{Lo, Li, L2) is the number of configurations if C Z"* n such that 
\K\ — n and Li < £i{K, r) < Li + 1 ioi i G {0, 1, 2}. This number can be estimated 
by the asymptotic value ^(^"1) Y[i=o-^i~^ when N 00. Therefore, taking (3 
large enough, we obtain that the first term of the RHS of (|93)) is bounded by 



Lj 2 
V V r7p-C"c2(3)|r| f 1±M^ ^ 



)^ - l|Im(z)| 



2+d 



With the same techniques, we find that the second and third terms of the RHS 
of ([M]) are respectively bounded by 



□ 



^\\lm{z)\J ^\\lm{z)\ 
which concludes the proof of the proposition. 
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